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Preface 


i 

During  a discussion  with  my  thesis  advisor  at  the  beginning  of  this 
investigation,  I first  mentioned  the  idea  of  utilizing  a rotating,  body 
fixed  coordinate  system  as  a possible  solution.  If  I had  gotten  that 
idea  at  any  other  time,  it  probably  would  have  died.  My  advisor  caught 
the  significance  of  that  particular  solution  method,  however,  and  he  was 
able  to  point  me  down  the  investigative  path  that  resulted  in  this  thesis. 

I don't  know  of  any  other  advisor  who  was  so  willing  to  spend  as  much  time 
discussing  results,  reviewing  drafts,  or  resolving  problems  with  his  students 
as  he  was.  I feel  quite  fortunate  to  have  had  Captain  James  E.  Negro  as 
my  thesis  advisor.  Thanks,  Jim. 

J.  G.  T. 
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Abstract 

A new  target  acceleration  model  capable  of  improving  the  performance  of 
an  airborne  aided  track  system  during  periods  of  tracker  sensor  failure  was 
investigated.  An  aided  track  target  prediction  algorithm  utilizing  range 
and  angle  rate  measurements  was  modified  so  that  the  target  state  estimates 
were  propagated  by  integrating  the  estimated  constant,  body  fixed  acceleration. 

The  modified  algorithm  performance  was  compared  to  the  original  by 
simulating  three  target  paths.  At  specified  times,  all  tracking  sensors 
were  turned  off  to  simulate  sensor  failure.  Each  algorithm  attempted  to 
keep  the  target  centered  on  the  tracker  boreslght  by  extrapolating  target 
position  estimates  without  processing  any  measurements.  roiS  pointing  errors 
were  computed  for  the  one  second  periods  following  each  failure. 

During  normal  sensor  operations,  the  performance  of  the  modified  algorithm 
was  equivalent  to  the  original.  When  the  sensors  were  failed  during  the 
tracking  of  targets  following  a constant,  body  fixed  acceleration,  the  modi- 
fied algorithm  generally  produced  smaller  RMS  errors.  (After  2.8  seconds 
of  filter  operating  time,  the  modified  algorithm  produced  one  third  the  RMS 
error  of  the  original.)  When  sensor  failure  occured  during  the  tracking 
of  targets  following  a rapidly  oscillating  path  or  traveling  in  a straight 
line  under  large  linear  accelerations,  both  algorithms  produced  equivalent 
RMS  errors.  The  modified  algorithm  demonstrated  the  same  sensitivity  to 
measurement  uncertainties  as  the  original. 
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IMPROVING  AIDED  TRACK  PERFORMANCE  DURING 
PERIODS  OF  TRACKER  SENSOR  FAILURE  BY  UTILIZING 
A TARGET  BODY  FIXED  COORDINATE  SYSTEM 


I . Introduction 


Background 

The  accuracy  requirements  that  exist  for  present  day  pointing  and 
tracking  systems  have  become  so  precise  that  conventional  feedback 
controllers  have  not  been  able  to  solely  satisfy  them.  One  system  of 
Interest,  used  to  track  airborne  vehicles,  uses  Kalman  filter  estimation 
techniques  to  improve  the  system  tracking  error  characteristics.  This 
system  produces  real  time  estimates  of  the  position,  velocity,  and  accel- 
eration of  a target.  These  estimates  are  then  processed  and  used  as 
additional  command  signals  to  the  tracker  glmbal  controls.  These  signals 
are  called  aiding  signals,  and  they  are  designed  to  minimize  the  dynamic 
lag  following  errors  that  occur  whenever  a target  is  tracked.  (Ref  1:25-38) 

Problem  Statement 

The  particular  tracking  system  described  in  Ref  1 is  designed  to 
operate  with  the  sensors  periodically  providing  discrete  target  Information 
(angle  rate  and  range)  to  the  aided  track  algorithm.  (This  condition  will 
be  refered  to  as  the  normal  mode.)  However,  the  target  tracking  sensors 
do  not  always  regiilarly  provide  information  because  they  cannot  alwavs 
detect  the  target.  Such  conditions  could  be  caused  by  momentary  sensor 
failure,  various  target  orientations,  jamming,  or  atmospheric  conditions. 
During  these  periods  of  sensor  blackout,  current  information  about  the 
target  is  no  longer  available  to  the  tracker  gimbal  controller.  However, 
the  system  can  utilize  the  last  information  that  it  processed  before 
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sensor  blackout  and  extrapolate  estimates  of  the  target  position,  velocity, 
and  acceleration  to  future  times.  (Such  a procedure  will  be  refered  to 
as  the  coast  mode . ) 

During  the  coast  mode,  commands  can  still  be  provided  to  the  gimbal 
controls  in  an  attempt  to  keep  the  target  aligned  with  the  boresight  of 
the  inoperative  sensors.  If  an  aided  track  algorithm  is  designed  that 
provides  a better  estimate  of  target  motion  during  the  coast  mode,  then 
the  system  will  keep  the  target  more  accurately  aligned  with  the  sensor 
boresight.  When  the  coast  mode  ends  and  the  sensors  operate  again,  the 
better  estimate  of  the  true  target  motion  will  result  in  a smaller  pointing 
error  and  a quicker  reacquisition  of  the  target.  This  thesis  develops 
and  investigates  the  performance  of  an  aided  track  algorithm  modification 
that  improves  the  performance  of  this  particular  tracking  system  during 
periods  of  sensor  blackout. 

Proposed  Solution 

The  performance  of  the  system  to  be  modified  during  this  investigation 
relies  heavily  on  the  simplified  dynamic  model  used  to  represent  the  target 
in  the  aided  track  algorithm.  This  model  represents  the  target  as  a 
point  mass  possessing  inertial  position,  velocity,  and  acceleration  vectors. 
These  three  vectors  completely  describe  the  dynamic  state  of  the  target. 

The  present  algorithm  model  also  assumes  that  the  inertial  acceleration 
of  the  target  remains  constant.  The  future  position  and  velocity  are 
estimated  by  simply  integrating  forward  that  constant  acceleration. 

Aerodynamic  vehicles,  however,  do  not  generally  possess  constant 
inertial  accelerations.  Rather,  the  accelerations  experienced  by  aircraft 
and  airborne  missiles  are  due,  for  the  most  part,  to  specific  lift  (as  a 
result  of  the  operation  of  lifting  and  control  surfaces)  and  specific 
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lift  (through  the  application  of  thrust  or  drag  producing  devices).  These 
accelerations  are  closely  associated  with  the  target  body  fixed  axes. 

The  algorithm  studied  in  this  thesis  utilizes  the  approximation  that  the 
non-gravitational  target  acceleration  remains  constant  with  respect  to 
a target  body  fixed  coordinate  system. 

Investigation  Method 

The  solution  to  the  problem  posed  by  this  thesis  was  investigated  by 
computer  simulation.  A computer  program,  described  in  Ref  2,  that  modeled 
the  time  response  of  a tracking  system  was  run  on  the  CDC  6600  computer 
located  at  Wright-Patterson  Air  Force  Base,  Ohio.  This  program  modeled  the 
physical  and  electronic  parameters  of  the  tracking  system  as  well  as  the 
time  behavior  of  the  target,  tracker  base,  and  disturbance  inputs.  (Ref  3:7) 
The  program  was  modified  during  the  thesis  investigation  to  incorporate 
different  target  trajectories,  to  modify  the  aided  track  algorithm  of  the 
present  system,  to  auid  a coast  mode  simulation,  to  simulate  range  and  point- 
ing error  uncertainties  and  gyro  drifts,  and  to  add  a data  output  section. 

Thesis  Organization 

This  thesis  is  divided  into  six  chapters  and  includes  two  appendices. 
Chapter  II  states  the  assumptions  made  during  the  study  to  simplify  cal- 
culations and  limit  the  scope  of  the  investigation.  Chapter  III  describes 
the  original  algorithm  and  then  mathematically  develops  the  modified  algorithm 
based  on  the  body  fixed  acceleration  assumption.  Chapter  IV  describes  how 
the  investigation  was  conducted  and  Chapter  V presents  tabular  and  graphical 
results  of  the  investigation.  Chapter  V also  discusses  those  results  in 
detail.  Conclusions  and  recommendations  for  future  study  are  listed  in 
Chapter  VI.  Appendix  A provides  a detailed  account  of  how  the  original 
computer  simulation  program  was  modified  during  the  investigation. 
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Appendix  B provides  a detailed  description  of  the  three  different  target 
paths  simulated  during  this  thesis  study. 
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II.  Scope  of  Investigation  and  Assumptions 


During  the  investigation  of  this  thesis,  a number  of  assumptions  were 
made  in  order  to  limit  the  scope  of  the  investigation.  The  tracker  base 
was  assumed  to  be  motionless  in  all  cases,  so  that,  for  analysis  purposes, 
the  system  dynamics  would  be  driven  solely  by  target  motion.  Also  during 
this  study,  the  effects  of  acceleration  due  to  gravity  were  neglected. 

These  first  two  assumptions  were  made  because  during  the  operation  of 
the  actual  system,  both  base  motion  and  gravity  can  be  con^juted  and  their 
effects  accounted  for. 

Coast  modes  one  second  in  duration  were  Investigated,  The  target 
was  assumed  to  be  a point  target  so  that  tracker  pointing  errors  caused 
by  the  system  locklng-on  to  different  parts  of  a target  would  not  be 
considered.  Also,  only  a complete,  simultaneous  loss  of  range,  angle  rate, 
and  pointing  error  sensors  was  considered.  Intermittent  or  partial  sensor 
failure  was  not  investigated.  It  was  also  assumed  that  any  sensor  failures 
would  be  detected  by  the  system,  so  that  erroneous  measurements  due  to  any 
failures  would  not  be  included  in  the  target  state  estimation  calculations. 

A real  target  can  change  acceleration,  within  certain  physical  limits, 
in  an  unpredictable  fashion.  If  this  occurs,  the  target  trajectory  will 
also  be  changed  unpredictably.  After  a sensor  blackout,  any  unmodeled 
acceleration  changes  to  a trajectory  even  perfectly  known  prior  to  black- 
out will  cause  pointing  errors.  Rather  than  model  the  possible  effects  of 
these  uncertainties,  it  was  assumed  that  any  state  estimation  errors  caused 
by  acceleration  changes  after  a sensor  blackout  would  equally  affect  the 
pointing  error  performance  of  the  present  and  proposed  algorithms.  There- 
fore, the  target  trajectories  that  were  produced  by  the  computer  simulation 
during  this  Investigation  were  deterministically  rather  than  randomly 
generated . 
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Finally,  the  sole  performance  criterion  for  comparing  algorithms  in 
this  study  was  root  mean  squared  (RMS)  tracker  pointing  error.  If  the 
RMS  pointing  error  produced  by  one  algorithm  during  a one  second  coast 
mode  period  was  smaller  than  that  produced  by  the  other  algorithm  for  the 
identical  period,  then  the  first  algorithm  was  assumed  to  be  better.  The 
instantaneous  error  present  at  the  end  of  a coast  mode  period  is  more 
important  than  RMS  error  because  of  the  major  effect  such  an  error  would 
have  on  the  reacquisition  dynamics  of  the  system.  However,  the  reacqtii- 
sitlon  time  of  the  system  will  not  occur  at  a precisely  known  time,  so  the 
RMS  error  computed  over  the  whole  coast  mode  period  will  give  a better, 
overall  measure  of  the  performance  for  a range  of  possible  reacquisition 
times. 
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III.  Algorithm  Modification 


Original  Algorithm 

Overall  System  Description.  To  place  the  following  discussion  of  the 
aided  track  algorithm  into  context,  a short  description  of  the  overall 
tracking  system  is  provided  in  this  subsection. 

The  angular  rate  command  signal  (azimuth  or  elevation)  that  is  provided 
to  the  respective  gimbal  actuator  is  the  sum  of  two  separate  signals.  One 
of  these  signals,  called  the  tracker  controller  signal,  is  produced  by  a 
proportional  plus  integral  feedback  scheme  that  uses  the  angular  error 
detected  by  a tracker/imager  as  the  feedback  error  signal.  (Ref  2:3-5) 

The  tracker  controller  provides  the  only  tracking  commands  until  the  aided 
track  algorithm  is  activated. 

The  other  part  of  the  angular  rate  command,  called  the  aiding  signal, 
is  produced  by  the  aided  track  algorithm  as  it  processes  range  and  angular 
rate  measurements.  The  aided  track  algorithm  is  described  in  greater 
detail  in  the  remainder  of  this  section.  A simple  functional  block 
diagram  describing  the  control  system  is  provided  in  Fig  1 . 

Propagation.  The  original  aided  track  algorithm  propagates  the 
estimates  of  the  target  state  vectors  from  one  update  time  to  the  next.  The 
dynamic  model  that  this  algorithm  employs  was  based  on  the  assumption  that 
the  target  is  a point  mass  that  has  associated  with  it  three  estimated 
state  vectors;  Inertial  acceleration,  inertial  velocity,  and  Inertial 
position.  The  estimated  inertial  acceleration  is  assumed  to  be  constant 
and  the  velocity  and  position  estimates  are  propagated  by  simply  integrating 
forward  that  constant  acceleration.  Mathematically  this  can  be 
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AIDED  TRACK  ALGORITHM 


Figure  1 . Simplified  Functional  Block  Diagram 
of  One  Channel  of  the  Control  System 
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expressed  as 


\(t  + AT) 

1 AT  ^ 

\(t) 

Vj^(t  + at) 

= 

0 1 at 

\(t  + AT) 

0 0 1 

(1) 


where  k = 1,  2,  3 (the  three  axes  of  an  inertial  Catresian  coordinate 
system),  t is  the  current  time,  AT  is  the  time  increment  over  which  the 
target  states  are  being  propagated,  and 
components  of  the  inertial  position,  velocity,  and  acceleration  vectors, 
respectively.  (Ref  1:32) 

The  current  estimates  of  position,  velocity,  and  acceleration  are 
converted  through  a coordinate  transformation  into  estimates  of  the  req.uired 
gimbal  azimuth  and  elevation  rates.  The  aided  track  algorithm  then  produces 
aiding  signals  that  are  added  to  the  tracker  controller  commanded  signals. 


, Vj^,  aJid  Aj^  are  the  k 


(Ref  li5) 

Updating.  At  each  update  time,  angle  rate  and  range  measurements  are 
taken.  These  measurements  are  transformed  into  a rectangular  set  of  (so 
called)  pseudo  measurements  which  consist  of  three  components  of  position 
and  three  components  of  velocity.  (Ref  1S30-31)  The  Kalmam  filter 
combines  the  pseudo  measurements  with  the  current  estimates  and  produces 
aji  updated  set  of  estimated  target  state  vectors.  (Ref  1:33)  These  updated 
estimates  are  then  propagated  by  the  aided  track  algorithm  to  the  next 
update  time.  The  fundamental  matrix  of  this  algorithm  models  the  acceler- 
ation of  the  target  as  a random  walk  phenomenon.  (Ref  1:3^) 

Coast  Mode.  The  original  computer  simulation  did  not  provide  for  a 
coast  mode  of  operation.  As  originally  programmed,  when  the  simulated 
target  could  no  longer  be  detected  by  the  tracker  sensors,  the  entire 
simulation  was  terminated.  (Ref  2:35”36)  The  program  was  modified, 
therefore,  so  that  this  important  phase  of  operation  could  be  investigated. 
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This  subsection  describes  the  performance  of  the  original  system  after 
it  was  modified  to  allow  the  simulation  of  the  coast  mode. 


Just  as  in  the  normal  mode,  the  estimated  state  vectors  are  propagated 
to  the  next  update  time.  However,  no  measurements  can  be  taJcen  at  that 
time  because  of  the  sensor  blackout,  and  the  Kalman  gain  equations  are 
not  processed  to  avoid  calculations  with  erroneous  data.  Instead,  updating 
does  not  take  place  and  the  state  estimates  simply  remain  unchanged.  This 
process  of  propagating  without  updating  is  continued  until  the  coast  mode 
ends  ajid  normal  sensor  measurements  are  again  available. 

During  a sensor  blackout  the  pointing  error  signals  obtained  from  the 
tracker/imager  are  unreliable.  At  the  beginning  of  the  coast  mode,  there- 
fore, the  pointing  error  signals  are  set  to  zero.  This  rate  memory  mode 
of  operation  holds  the  tracker  controller  signal  constant  at  a value  equal 
to  the  last  integrator  output  value.  The  aided  track  algorithm  continues 
to  provide  an  aiding  signal  based  on  the  current,  nonupdated  estimates 
of  the  target  state  vectors. 

Modified  Algorithm 

Modification  Rationale.  The  original  algorithm  assumed  that  the 
acceleration  of  am  aerodynamic  vehicle  (caused  by  the  forces  of  thrust, 
list,  and  drag)  is  inertlally  constant.  (This  analysis  neglects  gravity.) 
However,  in  the  aerodynamic  vehicles  that  will  be  tracked  by  this  system 
thrust  is  produced  by  engines  fixed  in  the  vehicles.  In  addition,  lift 
and  drag  are  defined  as  the  perpendicular  forces  resulting  from  aerodynamic 
effects  on  the  vehicles.  The  directions  of  these  aerodynamic  forces  are 
defined  with  respect  to  the  velocities  of  the  vehicles,  not  inertial  space. 
As  the  vehicles  maneuver,  the  orientations  of  thrust  and  these  aero- 
dynamic forces  also  change  with  respect  to  inertial  space. 
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These  forces  are  not  necessarily  constant  with  respect  to  a target 
fixed  coordinate  system,  but  they  are  more  closely  associated  with  such 
a system  than  with  an  inertially  fixed  system,  especially  if  the  target 
is  maneuvering.  Such  an  observation  about  the  forces  that  act  on  aircraft 
and  airborne  missiles  suggests  that,  to  improve  tracker  performance,  the 
original,  algorithm  could  be  modified.  The  resulting  modification  is  based 
on  the  assumption  that,  neglecting  the  acceleration  caused  by  gravity,  the 
magnitude  and  direction  of  the  acceleration  of  an  airborne  vehicle  are 
constauit  with  respect  to  the  vehicle  itself. 

Because  a pilot  (or  missile  controller)  can  change  the  throttle 
setting  or  angle  of  attack  of  a target  vehicle  at  any  time,  the  target  body 
fixed  acceleration  will  generally  not  be  constant.  However,  if  such  changes 
are  assumed  to  be  pertubations  in  a constant,  target  body  fixed  acceleration, 
then  an  estimate  of  that  acceleration  can  be  produced  and  propagated  by 
the  Kalman  filter.  The  resviltlng  target  state  estimates  can  be  trans- 
formed into  aiding  signals  and  added  to  the  tracker  controller  signals  as 
in  the  original  system. 

Because  this  model  represents  an  actual  aircraft  more  realistically 
than  the  original  model,  it  is  reasonable  to  expect  better  performance 
during  the  coast  mode  from  the  modified  algorithm  if  an  aircraft  is  tracked. 
The  remainder  of  this  thesis  develops  and  investigates  this  proposed 
modification. 

Analytic  Development.  To  tedce  advantage  of  the  knowledge  of  the 
dynamics  of  the  target  vehicles,  a target  model  with  a constant,  body 
fixed  acceleration  was  developed  to  be  used  in  the  propagation  portion  of 
the  aided  track  algorithm.  As  the  new  target  model  rotates,  the  estimated 
acceleration  vector  is  also  rotated  about  the  same  axis,  and  the  magnitude 
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of  the  acceleration  remains  constant.  (This  model  happens  to  represent 
perfectly  an  aircraft  performing  a level,  constant  rate,  constaint  airspeed 
turn.) 


Information  about  the  target  orientation  is  not  available  to  the 
algorithm  from  the  tracker  sensors.  Therefore,  the  longitudinal  axis  of 
the  target  model  was  defined  in  the  modified  algorithm  to  coincide  with 
the  estimated  velocity  vector.  A second  coordinate  direction  was  defined 
to  coincide  with  the  vector  cross  product  of  the  estimated  velocity  and 
estimated  acceleration  vectors.  The  third  body  fixed  coordinate  was  defined 
to  complete  a right  handed,  orthogonal  set. 

These  directions  maJce  up  the  U coordinate  system  which  is  defined  at 
any  time  t by  the  following  unit  vectors. 


V(t) 

(2) 

v(t)  X A(t) 
lv(t)XA(t)| 

(3) 

\(t)  = 1:^(4)  X s,(t) 

w 

where  the  superior  bar  represents  vector 
inertial  velocity,  A(t)  is  the 
estimated  inertial  acceleration, 

X is  the  cross  product  operator, 

Uy(t)  is  a unit  vector  in  the 
direction  of  V(t),  is  a 

unit  vector  in  the  direction  of 
V(t)  X A(t),  and  Uj^(t)  is  a unit 


quantities,  V(t)  is  the  estimated 


vector  that  completes  the  set. 
Fig  2 Illustrates  the  U system. 
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Figure  2.  U Coordinate  System 


Propasation.  The  purpose  of  the  propagation  section  of  this  algorithm 
is  to  find  the  estimated  future  values  of  A(t  + AT) , V(t  + AT),  and 
R(t  + aT)  given  the  current  best  estimates  of  A(t),  V(t),  and  R(t). 

The  following  subsection  of  this  thesis  develops  the  equations  to  find 
those  estimated  future  values. 

The  estimated  inertial  acceleration  and  velocity  vectors  of  the  target 
can  be  defined  at  any  time  t in  terms  of  the  rotating  U coordinate  system 
and  scalar  quantities  as 

A(t)  = Ay(t)^(t)  + Aj^(t)u^(t)  (5) 

V(t)  P V(t)i;y(t)  (6) 

These  vectors,  along  with  R(t),  the  estimated  position  vector,  are 
illustrated  in  Fig  3* 


point  target 


tracker 


V(t) 

v(t 


(U  coordinate  system) 


stationary 
ba 


Figure  3*  State  Vectors  Ebcpressed  in 
U Coordinate  System 


The  modified  algorithm  propagates  the  state  estimates  by  assuming 
that  the  components  of  acceleration  remain  constant  during  the  short  time 
period,  aT,  so  that  Ay(t)  = A^{t  + aT)  and  A^(t)  = A^(t  + aT).  The 
estimated  acceleration  and  velocity  vectors  for  time  t + aT  are  then 
found  to  be 


A(t  + at)  = Ay(t)uy(t  + at)  + Aj^(t)u^(t  + at) 
V(t  + at)  = V(t  +AT)u^(t  + at) 


(?) 

(8) 


The  scalar  V(t  + aT)  can  be  found  by  vector  calculus  as  follows. 

The  acceleration  vector  can  be  represented  by  applying  the  Theorem  of 
Coriolis  to  V(t). 


PjV(t)  = PyV(t)  + Wjy(t)  X V(t) 


(9) 


where  I represents  an  inertial  coordinate  system,  p represents  the  time 
derivative  operator  in  the  subscripted  coordinate  system,  and  W^y(t)  is  the 
angular  rate  of  the  U system  with  respect  to  the  I system.  Expressed  in 
the  U coordinate  system,  Wjy(t)  is 


V*) ' "lu  ^ ”iu 

V L N 


(10) 


If  equation  (9)  is  coordinatized  in  the  U system,  the  following  is  obtained; 


u 

V(t) 

u 

AL(t) 

= 

0 

+ 

0 

0 

lU 


u 

v(t) 

u 

(t) 

0 (t) 

0 

(11) 

(t) 

1 

1 

O 

M 

0 

where  the  superscripts  represent  the  reference  system  of  coordinatization, 
the  superior  dot  represents  the  scalar  time  derivative,  and  the  3 by  3 
matrix  is  the  skew  symmetric  form  of  Wjy(t).  (Ref  4t22) 
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From  matrix  equation  (ll),  the  following  scalar  equations  can  be 


obtained ; 


V(t)  = A^,(t) 

(12) 

= \t) 

(13) 

^lU  (^)^(^)  ~ ^ 

L 

(14) 

Since  (t)  does  not  appear  in  these  scalar  equations,  it  cannot  be  found 

lUy 

mathematically  from  them.  This  component  of  Wjy(t)  represents  target  roll. 
There  are  many  possible  methods  to  model  target  roll,  but  in  this  development 
target  roll  will  be  assumed  to  be  zero  so  that  (t)  = 0.  Also,  because 

lUy 

in  general  V(t)  / 0,  equation  (14)  results  in  (t)  = 0.  Therefore 

L 

“luf*-)  = “lu 

N 

Because  Ay(t)  has  been  assumed  to  be  constant,  simple  scalar  integration  of 
equation  (12)  from  time  = t to  time  t = t + aT  yields 

V(t  + at)  = Ay(t)AT  + V(t)  (16) 

At  this  point  in  the  derivation  of  the  new  propagation  equations,  it 
is  convenient  to  redefine  time  as  t = 0 at  the  beginning  of  each  propagation 
cycle.  Without  loss  of  generality,  the  inertial  coordinate  system  I, 
composed  of  unit  vectors  i,  j,  and  k,  can  also  be  defined  at  the  beginning 
of  each  propagation  cycle  so  that  the  U system  and  the  I system  coincide 
at  t = 0.  (uy(0)  = i,  Uj^(O)  = j,  ajid  Uj^(O)  = k)  The  unit  vectors  of  the 
rotating  U system  can  be  found  at  any  time  in  terms  of  the  non-rotating  I 
system  as 
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Uy(t)  = cos[0(t)]r  + sin[0(t)]j 

(17a) 

Uj^(t)  = -sin[9(t)Ji  + cos[0(t)]j 

(17t) 

(17c) 

where  9(t)  is  the  angle  through  which  the  U system  has  rotated  since 
t = 0. 

This  relationship  is  only  valid  for  0 £ t ^AT.  When  t = AT,  time 
is  redefined  as  t = 0 and  a new  I coordinate  system  is  defined  for  the 
next  propagation  cycle.  The  relationship  between  the  U and  I coordinate 
systems  is  illustrated  in  Fig  4. 


The  change  in  the  inertial  velocity  vector  that  occurs  from  time 
t = 0 to  time  t = AT  is  defined  as  the  vector  difference 

AV  = V(aT)  - V(0)  (18) 

Substitution  of  the  values  for  V(0)  and  V(aT)  from  equations  (6)  and  (8) 
into  equation  (18)  yields 

AV  = v(at)^(at)  - V(0)IIy(0)  (19) 
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Completing  the  transformation  defined  by  equation  (l?)  and  substituting 
for  V(aT)  from  equation  (l6)  gives  W in  terms  of  inertial  coordinates. 


AV  = [A^(0)aT  + V(0)][cos(^9(t)]i  + sin(^9(t)]j}  - V(0)i  (20) 

Substituting  into  equation  (I9)  and  solving  for  V(AT)uy(AT)  = V(aT)  yields 
V(aT)  = [Ay(0)AT  + V(0)]cos[9(t)]I  + [A^(0)aT  + V(0)]sin[9(t)]  j (21) 


9(aT)  can  be  fovind  by  integrating  equation  (13)- 


AT  at  k^{t) 

e(AT)  = S W (t)  dt  = S — dt  (22) 

0 0 V(t) 


V(t)  is  given  by  equation  (16)  and  Ay(t)  and  A^(t)  are  constants. 


Therefore 


, X AT  dt 

e(AT)  = A^(0)/  

0 Ay(0)t  + V(0) 


(23) 


Performing  the  indicated  integration  yields 

A,(0)  A^(0) 

0(aT)  = AT  + 1]  , A^(0)  / 0 (24) 

VO)  v(o)  ^ 

where  In  represents  the  natural  logarithm  operation. 

If  the  estimated  value  of  ky{0)  = 0,  the  value  for  9(aT)  becomes 
undefined  if  equation  (24)  is  used.  However,  if  Ay(0)  = 0,  then  V(t)  = V(0) 
and  the  estimated  magnitude  of  W^,  (t)  will  be  constant  during  time  AT. 
Subsequent  integration  of  equation  (23)  yields 

^(0) 

0(aT)  = ^ AT  , Av(0)  = 0 (25) 

v(0)  ’ 

which  is  to  be  used  in  place  of  equation  (24)  in  those  cases.  V(aT)  can 
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be  foimd  by  substituting  the  resulting  value  for  0(aT)  into  equation  (21). 

A(aT),  the  inertial  acceleration  vector,  can  be  found  more  directly, 
because  the  components  Ay(t)  and  are  assumed  to  be  constant.  Ihus, 

only  a coordinate  transformation  is  required.  Applying  the  coordinate 
transformation  defined  in  equation  (l?)  to  equation  (5)  and  gathering 
terms  resxilt  in 


A(aT) 


{Ay(0)cos[e(AT)]  - A^(O)sin[0(AT)]}l 

+ {Ay(O)sin[0(AT)]  + A^(O)cosC0(aT)]} j 


(26) 


R(aT)  can  be  foxond  by  integrating  the  components  of  V(t)  from  time 
t = 0 to  time  t = AT  and  adding  the  results  to  R(0).  However,  the  resulting 
integrands  contain  terms  of  the  form 


AT 

A.(0) 

Kio) 

I, 

= S 

sinC-^^ 

ln[ — — t + 

1]}  dt 

(27) 

1 

0 

V(0) 

AT 

A.(0) 

= S 

cos{ 

ln[  ^ t + 

-P 

1 — 1 

(28) 

2 

0 

V(0) 

which  have  not  been 

solved 

in  closed 

form. 

R(aT)  is  foimd  Instead  by  making  the  simplifying  assumption  that  during 
the  short  time  Af 


R(aT)  = R(0)  + V^aT  (29) 

where  is  a constant  velocity  vector  that  is  assumed  to  have  operated 
during  at. 

Because  V(0)  and  V(aT)  are  both  known,  the  assumption  is  made  that 

V,  . = V(0)  . ^ (30) 

where  aV  Is  given  by  equation  (20). 
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Substituting  into  equation  (29)  results  in 

R(aT)  = R(0)  + [V(0)  + ^]aT  (31) 

R(0)  is  known  in  terms  of  the  I system  as 

R(o)  = Ri(o)T  + Rj(o)j  + Rj^(o)k  (32) 

Substituting  into  equation  (3l)  for  R(0)  from  equation  (32)  and  f or  aV 
from  equation  (20),  then  gathering  terms  gives  the  desired  result: 

r(at)  = [r^(o)  + [A^,(0)aT  + V(0)]^^^^^^}^T]i 

+ (Rj(0)  + C[a,/o)at  + (33) 

+ CRk(0)]k 

If  during  ainy  propagation  cycle  the  estimated  value  of  either  Aj^(O) 
or  V(0)  is  zero,  then  the  definition  of  the  U coordinate  system,  equa^ions 
(2),  (3)»  and  (4),  cannot  be  applied.  In  those  cases,  the  modified  algorithm 
is  switched  back  to  the  original  algorithm  method  of  propagation.  (In  the 
first  case,  no  target  rotation  is  estimated,  and  in  the  second  case,  the 
target  is  estimated  to  be  motionless.  In  either  case,  the  original 
algorithm  is  expected  to  provide  satisfactory  propagation.) 

A summary  of  the  modified  algorithm  propagation  equations  studied  in 
this  thesis  is  listed  in  Fig  5* 

Updating.  The  equations  used  in  the  propagation  portion  of  the  modified 
algorithm  cannot  be  directly  incorporated  into  the  Kalman  gain  equations 
iised  in  the  update  portion  of  the  algorithm  because  they  are  non-linear. 

No  attempt  was  made  to  linearize  these  equations,  rather  the  original 
Kalman  gain  equations  from  the  unmodified  algorithm  were  retained  in  their 
original  form  so  that  their  practical  real-time  calculation  properties 
could  be  utilized. 
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Coast  Mode.  The  coast  mode  of  operation,  as  described  for  the  original 
algorithm,  was  also  unchanged  in  the  modified  algorithm.  The  same  procedure 
of  propagating  without  updating  was  followed  to  avoid  erroneous  measurements. 


A(t  + at)  = {Ay(t)cos[0(t  + at)]  - Aj^(t)sin[0(t  + AT)]}i 
+{Ay(t)sin[0(t  + at)]  + Aj^(t)cos[0(t  + AT)]}j 


V(t  + at)  = lk^it)AT:  + V(t)]cos[0(t  + AT)]i 

+ [AY(t)AT  + V(t)]cos[0(t  + AT)]3 

R(t  + at)  = [R^(t)  + 


(21) 


^ + CAy(t)AT  + v(t)]^°iCQ(y 
[Rj(t)  + iLi^{t)AT  ^ v(t)]gj£LQ(t  ^ at)]3^^i. 
+ C\(t)]k 


where 


or 


ArCt)  Ay(t) 

A,(t) 

0(t  + AT)  = , Ay(t)  = 0 


(24) 

(25) 


(aT,  the  estimated  values  of  Ay(t) , Aj^(t),  V(t),  R^(t),  Rj(t), 
Rj^(t),  and  directions  of  i,  j,  and  k are  known  at  the  beginning 
of  each  propagation  cycle.) 


Figure  5*  Summary  of  Propagation  Equations 
from  Time  t to  Time  t + aT 
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IV.  Investigation 

The  problem  posed  by  this  thesis  was  investigated  by  computer  simulation. 
The  CDC  6600  computer  located  at  Wright-Patterson  Air  Force  Base,  Ohio, 
was  used  to  execute  the  similLation  program  (described  in  Ref  2)  which 
modeled  the  time  response  characteristics  of  the  tracking  system.  This 
program  was  changed  to  incorporate  the  aided  track  algorithm  modifications 
that  were  developed  in  Chapter  III  of  this  thesis.  A coast  mode  simulation 
was  also  added  to  the  program  so  that  both  the  normal  and  coast  modes  of 
operation  could  be  studied.  Appendix  A provides  computer  listings  of  the 
modifications  to  the  original  computer  simulation  program. 

Method 

This  problem  was  investigated  by  simulating  the  tracking  of  three 
different,  representative  target  paths  and  comparing  the  performance  of  the 
original  system  with  that  of  the  modified  system.  The  coast  mode  initiation 
time  (designated  TV AN)  was  varied  during  each  computer  run,  so  that  any 
system  performance  sensitivity  to  shorter  or  longer  filter  operating  times 
could  be  found.  The  duration  of  the  coast  mode  was  1.0  seconds  for  all 
simulations.  The  effects  of  measurement  errors  were  Investigated  by 
simulating  range  and  pointing  error  uncertainties  and  gyro  drifts.  Diiring 
every  run  the  aided  track  algorithm  was  tvirned  on  at  t = 0.2  seconds.  The 
tracker  controller  alone  controlled  the  system  prior  to  that  time. 

RMS  azimuth  and  elevation  pointing  errors  were  computed  during  each 
one  second  coast  mode  period.  These  RMS  error  values  were  used  as  the 
sole  performance  criterion  during  this  study. 

Two  computer  runs  were  made  for  each  set  of  conditions.  One  run 
simulated  the  modified  system,*  the  other  simulated  the  unmodified  system. 
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The  two  runs  were  identical  except  for  the  insertion  of  the  algorithm 
modification  into  the  second  run.  Any  difference  in  performance  between 
the  two  simulated  runs  was  assumed  to  be  caused  by  the  algorithm  modification. 

Target  Paths 

The  target  paths  were  designated  Path  A,  Path  B,  and  Path  C.  Path  A 
simulated  a constant  speed,  circulair  target  motion;  Path  B simulated  a 
constant  speed,  sinusoidal  target  motion;  and  Path  C simulated  a straight 
line  target  motion  under  the  influence  of  a high  g,  one  dimensional,  inertial 
acceleration  followed  by  a similar  deceleration.  These  paths  are  illustrated 
in  Fig  6 £ind  are  described  in  detail  in  Appendix  B. 

Error  Simulation 

The  system  being  studied  obtains  measurements  from  three  sovirces.  The 
angular  rates  that  are  used  as  commands  to  the  gimbal  actuators  are  used 
as  the  system  inertial  boresight  angular  rate  measijrements , range  information 
is  obtained  from  a ranging  device,  and  tracker  boresight  error  is  obtained 
from  an  infra-red  tracker/imager  device.  (Ref  2:3)  These  real  devices 
are  affected  by  errors  that  the  original  simulation  program  did  not  model. 

In  order  to  investigate  the  sensitivity  of  the  modified  system  to  measurement 
errors  or  gyro  drifts,  the  computer  program  was  modified  to  model  those 
effects. 

Random  Measurement  Errors.  The  original  computer  program  provided 
a random  number  generator  that  was  utilized  to  model  the  random  measurement 
errors  that  occur  in  the  real  ranging  and  tracker/ imager  devices.  Three 
Independent  sequences  of  white,  Gai^ssian  random  numbers  with  zero  means 
and  specified  standard  deviations  were  generated  and  added  to  the  modeled 
measurements  of  elevation  and  azimuth  pointing  and  range  errors,  respectively. 
The  standard  deviations  were  set  at  5*0  meters  for  the  range  error  and 
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PATH  B 


constant  speed  = 600.42  m/sec 
(approx.) 


PATH  C 


V/.'i  _ r-300  m/secp,  t<2  sec. 
i'  ' C 40  m/sec  , ti2  sec. 
(speed  varies) 


Figvtre  6.  Illustration  of  Three  Investigated 
Target  Paths 
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10.0  mlcroradians  for  the  two  angular  measurement  errors.  It  was  assumed 
that  all  other  range  or  boresight  error  measurement  uncertainties  were 
small  when  compared  to  these  random  measurement  errors. 

To  provide  a common  basis  for  comparing  results,  the  same  set  of 
random  sequences  was  used  for  every  simulation  run.  To  insure  that  the 
choice  of  one  particular  set  of  random  sequences  did  not  influence  the 
final  results,  several  simulations  runs  were  made  with  no  uncertainty 
inputs.  By  comparing  the  results  of  simulations  with  uncertainty  inputs 
to  those  without,  it  could  be  determined  whether  this  specific  set  of 
uncertaintly  sequences  or  the  basic  difference  between  the  two  algorithms 
caused  the  results  reported  later  in  this  study. 

Gyro  Drift  Errors.  It  was  assumed  that  the  only  significant  contribution 
to  the  angular  rate  measurement  errors  would  be  caused  by  gyro  bias  drifts. 

To  simulate  these  errors,  a constant  value  of  0.1  degrees/hour  (approx- 
imately 0.48  mlcroradians/second)  was  added  to  the  modeled  elevation  and 
azimuth  inertial  rate  values  produced  by  the  section  of  the  computer  program 
that  simulated  the  gyro  dynamics  of  the  system. 
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V.  Results  and  Discussion 


Resvilts 

In  this  section  of  the  report  the  results  of  this  investigation  are 
presented  in  tabular  and  graphical  form.  In  each  of  the  following  tables 
and  graphs,  TVAN,  the  coast  mode  initiation  time,  is  the  varying  parameter. 
Separate  tables  and  graphs  are  presented  for  each  of  the  three  represent- 
ative target  paths.  Those  tables  and  graphs  are  grouped  according  to 
target  path. 

Tables.  Three  tables  of  results  are  presented.  Each  table  compares 
the  RMS  azimuth  and  elevation  pointing  errors  of  the  unmodified  system 
with  the  corresponding  errors  of  the  modified  system  during  six  different 
one  second  coast  mode  periods.  Each  coast  mode  period  was  initiated  at 
a different  TVAN,  starting  at  t = 0.5  seconds  and  increasing  by  O.5  second 
increments  until  t = 3.0  seconds.  The  RIB  errors  were  computed  by  the 
computer  program  for  the  period  t = TVAN  until  t = TVAN  +1.0  second  for 
each  run. 

Graphs . Six  graphs  of  results  are  presented.  Alternating  graphs  plot 
azimuth  or  elevation  pointing  error  versus  time  for  both  the  modified  and 
unmodified  systems,  one  for  each  target  path.  The  results  of  varying  TVAN 
for  each  target  path  are  shown  by  plotting  all  the  runs  for  the  various 
TVAN  times  on  the  same  graph.  Each  plot  ends  1.0  second  after  the  respective 
TVAN  time.  Four  additional  graphs  are  presented  in  the  Discussion  section 
of  this  chapter  which  illustrate  the  effects  of  the  measurement  uncertainties 
on  the  system  performance  and  which  Illustrate  the  basic  difference 
between  the  two  propagation  methods. 
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Table  I 


Coast  Mode  RMS  Angular  Pointing  Error  Comparison,  Path  A 


TVAN 

Chauinel 

Unmodified  Algorithm 

Modified  Algorithm 

(Seconds) 

RMS  Error  (Milliradians) 

RMS  Error  (Milliradians) 

0.5 

Azimuth 

4 .4092 

4.4094 

Elevation 

1.1522 

1.1555 

1.0 

Azimuth 

3.9983 

4.0010 

Elevation 

0.1036 

0.1459 

1.5 

Azimuth 

1.3632 

1.3590 

Elevation 

0.3165 

0.1788 

o n 

Azimuth 

0.4103 

0.2847 

Elevation 

0.1588 

0.0635 

2.5 

Azimuth 

0.9088 

0.4113 

Elevation 

0.8758 

0.3045 

3.0 

Azimuth 

1.0272 

0.3492 

Elevation 

1.6985 

0.6857 
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RZinUTH  POINTING  ERROR  (NRRDIflNS) 


PRTH  fl 


TIME  (SECONDS) 


Figure  ?•  Azimuth  Pointing  Error  versus  Time, 
Various  TVAN  Times,  Path  A 
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ELEVATION  POINTINO  ERROR  (MRROIANS) 
rlZ.OO  -8.00  -4.00  0.00  4.00  8.00  18.00 


PATH  fl 


TIME  (SECONDS) 


Figure  8.  Elevation  Pointing  Error  versus  Time, 
Various  TV AN  Times,  Path  A 
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Table  II 


Coast  Mode  RMS  Angular  Pointing  Error  Comparison,  Path  B 


TVAN 

Channel 

Unmodified  Algorithm 

Modified  Algorithm 

(Seconds) 

RMS  Error  (Milliradians) 

RMS  Error  (Milliradians) 

0.5 

Azimuth 

3.8033 

3. 8033 

Elevation 

0.1516 

0.1527 

1 0 

Azimuth 

5.3652 

5.3647 

Elevation 

0.0579 

0.1265 

1.5 

Azimuth 

6.4504 

6.4449 

Elevation 

0.0184 

0.0419 

Azimuth 

10.2466 

10.2210 

Elevation 

0.0457 

0.4373 

2.5 

Azimuth 

11.6687 

11.6216 

Elevation 

0.1487 

0.0419 

3.0 

Azimuth 

11.5772 

11.6858 

Elevation 

0.0451 

0.1642 
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Figure  9«  Azimuth  Pointing  Error  versus  Time, 
Various  TV AN  Times,  Path  B 
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ELEVATION  POINTING  ERROR  (NRADIANS) 

,-30.00  -20.00  -10.00  0.00  10.00  20.00  30.00 


.00  1.00  2.00  3.00  4.00 


TIME  (SECONDS) 


Figure  10.  Elevation  Pointing  Error  versus  Time, 
Various  TV AN  Times,  Path  B 
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Table  III 


Coast  Mode  RMS  Angiilar  Pointing  Error  Comparison,  Path  C 


TVAN 

Channel 

Unmodified  Algorithm 

Modified  Algorithm 

(Seconds) 

RMS  Error  (Milliradians) 

RMS  Error  (Milliradians) 

0.5 

Azimuth 

3.1699 

3.1700 

Elevation 

3.1997 

3.2002 

1 n 

Azimuth 

3.7^3 

3.7392 

1 tU 

Elevation 

3.5835 

3.4873 

Azimuth 

1.8031 

1.7459 

1.5 

Elevation 

1.8103 

1.6521 

2.0* 

Azimuth 

5.8712 

5.8629 

Elevation 

5.3699 

5.2887 

2.5 

Azimuth 

IO.3085 

10.4500 

Elevation 

9.1538 

9. 1089 

3.0 

Azimuth 

9.4819 

9.0621 

Elevation 

7.8576 

7.6997 

♦Deceleration  Time 
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Flgiire  11.  Azimuth  Pointing  Error  versus  Time, 
Various  TVAN  Times,  Path  C 
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ELEVATION  POINTING  ERROR  (MRADIANS) 
rSO.OQ  -20.00  -10.00  0.00  10.00  20.00  30.00 


.00  1.00  2.00  3.00  4.00 
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Figure  12.  Elevation  Pointing  Error  versus  Time, 
Various  TVAN  Times,  Path  C 
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Discussion 


This  thesis  investigated  improving  the  performance  of  an  aided  pointing 
and  tracking  system  during  periods  of  sensor  blackout.  It  is  significant, 
however,  than  during  the  normal  modes  of  operation,  prior  to  t = TVAN  in 
all  the  computer  simulations  performed,  the  pointing  error  performance  of 
the  modified  system  was  virt^lly  identical  to  that  of  the  original  system. 
This  result  was  obtained  because  during  the  short  propagation  time  period 
of  0.1  seconds,  the  two  different  aiding  algorithms  produced  only  slightly 
different  target  state  estimates.  These  slightly  different  estimates 
produced  only  slightly  different  aiding  signals  which,  when  added  to  the 
tracker  controller  signals,  produced  almost  identical  angular  rate  commands. 
These  commands  subsequently  resulted  in  the  nearly  identical  pointing  errors. 

The  difference  in  performance  between  the  two  systems  depended  on  the 
amount  that  the  estimated  acceleration  vector  was  rotated  in  the  modified 
algorithm.  That  amount  of  rotation,  0(aT),  was  small  because  AT  was  only 
0.1  seconds. 

If  extremely  large  magnitudes  of  (t)  occiired  during  the  period  AT, 

-LUn 

large  differences  between  the  performances  of  the  two  algorithms  would  occur. 

However,  (t)  is  limited  in  magnitude  by  the  aerodynamic  and  structural 

limits  of  the  airborne  target  vehicles.  Such  limited  values  of  (t), 

N 

integrated  over  a AT  of  0.1  seconds,  are  expected  to  continue  to  produce 
small  values  of  0(aT). 

Significant  differences  did  occur  between  the  two  systems  during  the 
coast  modes  of  operation,  after  t = TVAN.  The  following  subsections 
discuss  these  differences  at  they  occured  during  each  target  path  simulation. 

Path  A Discussion.  The  modified  system  generally  performed  much  better 
than  the  original  when  Path  A was  simulated.  This  result  was  to  be  expected, 
because  the  target  motion  that  was  simulated  by  Path  A was  exactly  the  type 
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of  motion  that  was  modeled  in  the  propagation  section  of  the  modified 
algorithm.  The  modification  did  not  produce  exact  target  tracking  during 
the  coast  mode  periods,  however,  because  the  uncertainties  in  the  measure- 
ments prevented  the  zero  error  target  state  estimates  from  being  calculated. 

Table  I illustrates  that  the  degree  of  improvement  shown  by  the  modified 
system  is  sensitive  to  the  filter  operating  time.  After  1.3  seconds  or 
more  of  filter  operating  time  (TVAN  =1.5  seconds  or  higher),  the  modifi- 
cation consistantly  produced  smaller  RMS  pointing  errors  in  both  channels. 
After  2.8  seconds  of  filter  operating  time  (TVAN  = 3,0  seconds)  the  modifi- 
cation produced  RMS  pointing  errors  that  were  approximately  one  third  the 
size  of  those  produced  by  the  original  algorithm. 

Prior  to  TVAN  =1.5  seconds,  the  modified  system  produced  RMS  pointing 
errors  that  were  the  same  to  slightly  larger  than  those  that  the  original 
system  produced.  At  the  earlier  times,  neither  aiding  algorithm  has  had 
enough  operating  time  to  calculate  accvirate  target  state  estimates.  The 
modified  algorithm  propagation  of  these  inaccurate  estimates  created  from 
equivalent  to  slightly  larger  errors  than  the  original  algorithm  propagation 
of  those  same  Inaccurate  estimates.  The  largest  amount  by  which  the  original 
system  performed  better  than  the  modified  system  occured  at  TVAN  =1.0 
seconds,  when  the  unmodified  elevation  channel  showed  an  RMS  error  appox- 
imately  0.0423  milliradians  smaller  thauri  the  modified  elevation  channel. 

At  all  other  times  prior  to  TVAN  =1.5  seconds,  the  performance  of  the  two 
systems  were  equivalent. 

The  reason  for  these  results  is  graphically  illustrated  in  Fig  13. 

This  flgiire  contains  two  graphs  that  plot  the  actual  and  estimated  values 
for  the  inertial  acceleration  versus  time  for  each  algorithm  as  Path  A 
is  tracked,  TVAN  = 3*0  seconds  for  these  plots  and  measurement  uncertainties 
are  Included.  Prior  to  t = TVAN  during  the  normal  mode,  the  estimated 
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values  for  X^Ct)  computed  by  each  algorithm  are  virtually  identical.  At 
t = TVAN  = 3«0  seconds,  however,  when  the  coast  mode  is  initiated,  the 
difference  between  the  two  algorithms  becomes  apparent. 

Because  the  unmodified  algorithm  models  a constant  inertial  accel- 
eration, the  estimated  value  of  X2(t)  remains  constant  after  t = TVAN. 

The  modified  algorithm,  on  the  other  hauid,  models  a rotating  inertial 
acceleration  vector,  so  after  t = TVAN  the  estimated  value  of  X2(t) 

(one  component  of  that  rotating  acceleration  vector)  propagates  as  a chang- 
ing acceleration.  In  the  case  of  Path  A,  this  changing  estimated  acceleration 
more  closely  resembles  the  actual  acceleration  and  the  RMS  errors  are  smaller 
as  a result.  The  characteristics  of  the  acceleration  estimates  illustrated 
in  Fig  13  are  typical  of  the  X^  and  X2  acceleration  estimates  computed 
during  all  the  simulations  of  Path  A. 

Fig  14  contains  two  graphs  that  plot  the  actual  and  est?,:iated  values 
for  the  X2  inertial  acceleration  versus  time  for  each  algorithm  as  Path  A 
is  tracked.  TVAN  = 3.0  seconds,  but  unlike  Fig  13,  no  measurement  uncertain- 
ties are  included.  Comparison  of  Fig  14  to  Fig  13  illustrates  that  the 
difference  in  performance  foimd  between  the  original  and  modified  algorithms 
is  caused  by  the  different  propagation  methods  utilized  by  each.  The  same 
result  is  obtained  from  the  perfect  measurement  case  as  from  the  uncertain 
measurement  case  except  that  the  acceleration  estimates  are  slightly  less 
accurate  in  the  uncertainty  case.  Because  the  results  of  these  two  cases 
are  similar  and  because  any  other  set  of  uncertainty  sequences  with  the  same 
zero  means  and  standard  deviations  is  expected  to  be  very  similar  to  this 
set,  it  can  be  concluded  that  only  one  set  of  uncertainties  need  be  inves- 
tigated in  this  study  to  produce  a general  result.  Fig  14  is  typical  of 
the  perfect  measurement  simulations  performed  in  this  study. 
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X2  flCCEL  (M/S«*2) 

r60.00  -40.00  -20.00  0.00 


Figure  13-  Actual  and  Estimated  Magnitudes  of  X^(t), 
Path  A,  TVAN  = 3*0  Seconds,  for  Unmodified 
and  Modified  Algorithms,  Measurement  Uncertainties  Included 
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X2  flCCEL  (M/S«w2) 

-60.00  -40.00  -20.00  0.00 


unmodified  algorithm 

D Actual  Acceleration 
^ Estimated  Acceleration 


0.00  1.00  2.00  3.00 

TIME  (SEC) 


4.00 


MODIFIED  ALGORITHM 
O Actual  Acceleration 
^ Estimated  Acceleration 


0.00  1.00  2.00  3.00 

TIME  (SEC) 


4.00 


Figure  14.  Actual  and  Estimated  Magnitudes  of  X„(t), 

Path  A,  TVAN  =3*0  Seconds,  for  Unmodified 
and  Modified  Algorithms,  Measurement  Uncertainties  Not  Included 
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Path  B Discussion.  The  performance  of  the  modified  and  unmodified 
systems  during  the  simulation  of  Path  B were  equivalent.  Observation  of 
the  azimuth  errors  listed  in  Table  II  shows  that  neither  method  was  clearly 
superior.  At  different  times  each  showed  slightly  better  performance  than 
the  other,  but  no  trends  indicating  one  aiding  algorithm  better  overall  are 
found  in  the  results. 

The  rapidly  changing  target  states  of  Path  B could  not  accurately  be 
calculated  by  the  Kalman  filter  because  the  filter  bandwidth  was  too  narrow 
to  respond  satisfactorilly  to  the  O.5  Hertz  target  oscillations.  Consequently, 
the  state  estimates  used  to  compute  the  aiding  signals  were  grossly  inaccurate 
in  both  algorithms. 

When  a coast  mode  period  would  begin  at  t = TVAN,  the  two  algorithms 
propagated  those  erroneous  estimates  into  the  future  by  two  different  methods. 
The  state  estimates  obtained  from  each  algorithm  after  this  process  were 
different  from  each  other,  but  they  were  both  even  more  different  from  the 
true  target  states.  Also,  the  target  motion  encountered  during  the  simulation 
of  Path  E was  not  modeled  properly  in  either  algorithm  because  the  acceleration 
was  not  constant  in  either  inertial  or  body  fixed  coordinates.  Because  each 
algorithm  started  the  coast  mode  with  equally  poor  estimates  of  the  true 
target  states,  because  neither  algorithm  properly  modeled  the  target  dynamics, 
and  because  the  two  estimated  target  states  produced  by  the  different  algorithms 
were  much  more  closely  equal  to  each  other  than  to  the  true  target  state, 
the  resulting  RMS  azimuth  errors  produced  by  both  systems  were  nearly  equal. 

The  RMS  azimuth  errors  Increased  as  TVAN  increased.  This  occured  because 
as  the  target  moved  closer  to  the  tracker,  the  required  tracker  angular 
deflections  and  rates  needed  to  properly  track  the  oscillating  target 
increased . 
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The  RMS  elevation  errors  listed  in  Table  II  illustrate  the  measurement 
uncertainty  effects  on  the  operation  of  the  two  aiding  algorithms.  No 
target  motion  took  place  out  of  the  inertial  plane,  therefore  all  the 

elevation  errors  listed  were  caused  strictly  by  the  measurement  uncertainties, 
not  target  dynamics  or  system  dynamic  lag.  Those  uncertainties  created 
erroneous  target  motion  estimates  that  propagated  out  of  the  X^X2  plane, 
causing  elevation  rate  commands  which,  in  turn,  caused  elevation  channel 
errors.  In  the  elevation  channel,  neither  algorithm  showed  more  sensitivity 
to  the  measurement  uncertainties;  each  algorithm  showed  better  elevation 
channel  performance  at  different  TVAN  times. 

Path  G Discussion.  The  performance  of  the  two  systems  during  the 
simulation  of  Path  C were  equivalent.  Table  III  illustrates  that  at  various 
TVAN  times,  each  system  produced  RMS  pointing  errors  that  were  slightly 
smaller  than  those  produced  by  the  other.  However,  neither  system  showed 
an  overall  superior  performance. 

This  equivalent  performance  is  not  surprising,  because  the  target  model 
that  was  incorporated  into  the  original  algorithm  perfectly  modeled  the 
actiial  target  motion  (disregarding  the  instantaneous  acceleration  change  at 
t = 2.0  seconds).  Also,  the  modified  algorithm  was  designed  to  switch  to 
the  original  propagation  method  in  cases,  such  as  Path  C,  when  A^('t)  = 0. 

If  the  algorithms  had  produced  perfect  state  estimates,  the  modified 
algorithm  would  have  calculated  = 0,  and  the  original  method  portion 

of  the  modified  algorithm  would  have  been  switched  on.  However,  because 
of  the  measurement  uncertainties,  an  estimated  value  of  ^ 0 was 

calculated  and  the  modified  algorithm  propagated  the  current  state  estimates 
differently  than  the  original  algorithm. 

These  calculated  estimated  values  of  were  small,  though,  so  that 

the  estimates  resulting  from  the  modified  algorithm  propagation  were  close 
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to  those  resulting  from  the  original  propagation.  Thus,  the  aiding 
signals  and  subsequent  RMS  pointing  errors  were  close  to  the  same  values. 

The  smallest  RMS  pointing  errors  produced  by  both  systems  occured  during 
the  coast  mode  period  that  started  at  TVAN  = 1.5  seconds.  This  result  illus- 
trates that  both  algorithms  were  producing  increasingly  better  estimates  of 
the  -X^  acceleration  as  the  operating  times  of  the  filters  increased.  If 
the  acceleration  of  the  target  had  remained  constant,  it  is  expected  that 
both  algorithms  would  continue  to  produce  increasingly  better,  nearly  equiv- 
alent, state  estimates,  and  the  values  of  subsequent  one  second  coast  mode 
period  RMS  errors  would  show  even  better  performance. 
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VI . Conclusions  and  Recommendations 


This  chapter  lists  the  conclusions  formulated  during  this  Investigation 
and  recommends  areas  for  further  study. 

Conclusions 

1.  The  coast  mode  performance  of  the  particular  aided  pointing  and 
tracking  system  studied  In  this  thesis  can  be  Improved,  under  certain  con- 
ditions, by  making  a simple  software  modification  to  the  aided  track  algorithm. 
This  modification  Incorporates  a body  fixed,  rotating  coordinate  system  Into 
the  propagation  pozrtlon  of  the  aided  track  algorithm.  The  coast  mode  per- 
formance of  this  particular  modification  does  not  show  Improvement  under 

all  conditions,  but  the  resulting  performance  Is  always  at  least  equivalent 
to  the  original  system  performance. 

2.  The  software  modification  suggested  by  this  thesis  can  be  Incor- 
porated Into  the  system  without  degrading  normal  mode  performance  of  the 
system,. 

3.  The  siiggested  modification  has  the  same  sensitivity  to  measurement 
uncertainties  as  the  original  system. 
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Recommendations 


1.  A follow-on  study  to  this  investigation  should  be  performed  that 
utilizes  a more  sophisticated  propagation  model.  This  new  model  should  not 
be  constrained  merely  to  rotate  a constant  magnitude,  inertial  acceleration 
vector  about  a body  fixed  axis  as  the  model  in  this  investigation.  Rather 
the  new  model  should  allow  the  effects  of  a changing  body  fixed  acceleration 
(changed  by  target  roll  or  changes  in  thrust  or  angle  of  attack)  to  be 
propagated . 

2.  A follow-on  study  to  this  investigation  should  be  performed  that 
checks  the  validity  of  the  assumption  made  in  this  thesis  that  the  effects 
of  gravity  can  be  detected  and  compensated  for.  This  suggested  study  should 
proceed  by  investigating  the  errors  that  would  create  uncertainties  in  the 
estimates  of  the  true  magnitude  and  direction  of  the  acceleration  due  to 
gravity.  The  performance  of  the  pointing  and  tracking  system  as  it  processes 
these  erroneous  gravity  estimates  could  then  be  investigated. 

3*  This  Investigation  assumed  one  set  of  measurement  error  values. 

The  relative  sensitivity  of  the  aided  pointing  and  tracking  system  to 
different  set's  of  measurement  errors  was  not  studied.  A study  should  be 
performed  that  investigates  changes  in  performance  caused  by  changing  the 
quality  (and  thus  the  accuracy)  of  the  sensing  devices  used  on  the  system. 

4.  This  thesis  did  not  Investigate  the  effects  on  the  performance  of 
the  pointing  and  tracking  system  that  adding  or  removing  different  sensing 
devices  would  cause.  A follow-on  study  should  be  performed  that  studies 
different  measurement  devices  and  their  effects  on  the  system. 

5.  A follow-on  study  to  this  investigation  should  be  performed  that 
investigates  the  effects  of  individual  sensor  blackouts  on  the  system 
pointing  error  performance.  This  thesis  studied  only  the  effects  of 
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complete,  simultaneous  sensor  losses.  In  actual  practice,  however,  sensors 
may  fail  individually  so  that  some  valid  sensor  measurements  would  still 
be  available.  The  performance  of  the  pointing  and  tracking  system  under 
those  conditions  shoiild  be  studied. 
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Appendix  A 


Computer  Program  Modification  Listings 


The  modifications  to  the  original  computer  program  were  incorporated 
into  the  program  by  utilizing  the  CDC  UPDATE  utility  program  available 
on  the  GDC  6600  computer  located  at  Wright-Patterson  Air  Force  Base,  Ohio. 
The  following  listings  were  divided  into  functional  groupings  and  are 
presented  here  to  illustrate  exactly  how  the  original  program  (described 
in  Ref  2)  was  modified. 


Basic  Algorithm  Modification 

The  results  of  Chapter  III  of  this  thesis  were  coded  and  inserted  into 


the  original  program  as  illustrated  below. 

♦IMSERT  MDT9K.23 

OIMF’^SION  t/OXftO(l),UN(3),UL(3),UV(3) 

♦INSERT  (\nTRK.238 

CC:  COMPUTE  aNOLE  BETWEEN  PRESENT  EST  VEL  BNO  ACCEL  VECTORS 

00TAV=C(2,1)*C(3,1)*C(2,2)  'C(3,2)+C(2,  3>*C(3,3) 
ANA30=SQRT(C{3, 1 ) ♦♦E+C ( 3 ,2 ) ♦ *2 ♦C ( 3, 3) **2) 

IF(A3S(AMa30> .LT.l.OE-3) GO  TO  599 
VMAGOsSORTCCCB, 1 ) ♦* 2+C (2 ,2 ) *♦ 2 ♦C ( 2, 3 ) ♦ *2) 

IF(  AOS  (VNAGO)  .l.OE-3)  GO  TO  599 
CALPHA  = D0T  AV<  < AMAGO*VMAGO) 

ALPHa=ACOS (CALPHA) 

SALPHA=SIN(ALPMfl) 

IF (A9S  (SALPHA) . LT . 1. OE-3) 3 1 TO  5^9 
CCS  COMPUTE  UNIT  VECTOR  U-3U3-V  (UV)  IN  INERTIAL  COORDS 

00  5o0  1=1*3 

560  UV(I)=C(2,I) /VMAGO 

cc:  COMPUTE  Jn:t  vector  U-IUO-n  CJN)  in  INERTIAL  COORDS 

V0XA0(1)=C(2,2)‘C(3,3)-C  (2,3)*C(3,2) 

V0XA0( 2) =C (2,3) *0 ( 3, 1) -C (2 ,1) ♦C (3,3) 

V0XA0(  3)  =C(2, 1)  *0(3,2)  -C  (2  ,2)  ♦C(3,l) 
VOXaOM=SQRT(VOXao(l)**2+VO XAO(2)**2*VOXAO(3)**2) 

DO  573  1=1,3 

570  UN(I)=V0XA0(I)/V0XA0M 

CCC  COMPUTE  UNIT  VECTOR  U-SU9-L  (UL)  IN  INERTIAL  COORDS 

UL(1)  = UN(2)*UV(3)-UN(3)''  UV  (2) 

UL(2)=IJN(3)*'JV(1)-UN(1)*UV  (3) 

UL(3)=UN(1) *UV( 2)-UN(2)*UV (1) 


4? 


cc:  CO'^P'JTt  :0'^>'’0NENTS  OF  aCOEL  IN  U C lOPO  SYSTEM 

ftVM=A'iAr,o*CaLO-'  ft 
ALM  = AM  AGO^SftLP-^ft 
CCO  CO'^PUTE  TH-TA 

ir(ftPS (A7H) .LT.l .OE-6)GO  TO  53G 
ARGMT=  AVM*OTT/\/MftGO  + l 

C (CHECK  FOP,  N''GATI\/E  OP  TFRO  aRGUM-NT  qt  NATURAL  LOG) 

IPCARG^^T.LE.O'O  TO  53C 
THETA  = AL  i^ftLOGf ARGHT) /AVM 
GO  TO  GOG 

53D  THETA= ALM*DTT/7«ftGO 
590  CTH=COS(THETA) 

STH=SIM(THETft) 

CCC  COMPUTE  N-n  ACCEL  ESTIMATE 

DO  60  0 1 = 1, ■' 

so:  EX(  3, 1 ) = (AVM*r‘^w  -ALM*ST«)*  UV  (I ) ♦ ( A VM*  3 TH+AL  M^CT  H ) •UL  ( I ) 

CCS  COMPUTE  NE^J  \/EL  ESTIMATE 

VMftGN=  AVM*0TT  + 7*^AG0 
DO  610  1=1,3 

610  EX(2,I)=VMAGN»CTH*UV(I)+VMftGM*STH*UL(I  ) 

CCS  COMPUTE  NEf)  POS  ESTIMATE 

RI=  (VMAGO/2+  (A\/M*OTT+VM  AGO)  *CTH/2)  *011 
RJ=  ( (AVM*OTT+VMAGO)*STH /2)*OTT 
DO  650  1=1,3 

650  EX  (1,I)=RI»UV(I) +RJ*UL(I)  TCd,!) 

DO  860  1=1,3 
DO  660  J=l,3 

IVINC=3.* ( J-1 +(1-1.) ♦250 
X(IVINC) =EX (I, J) 

J\/INC=IVIN0  + 13 
660  X(JVINC)=6(I, J) 

GO  TO  210 
599  CONTINUE  ' 
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Path  Simxilatlon 


The  following  statements  simulated  each  target  path  by  specifying 
the  time  history  of  the  inertial  acceleration  to  be  applied  to  the  target. 
The  computer  simulation  integrated  these  accelerations  to  compute  the  target 
. jlocity  and  position  for  later  use  in  the  simulation.  Following  each  set 
of  program  statements  is  presented  the  initial  condition  matrix  for  that 
particular  path.  The  computer  variable  names  given  to  the  set  of 

inertial  target  states  is  given  by 


x^(t) 

X2(t) 

X3(t) 

X(09) 

X(12) 

X(15) 

x(t)  = 

Xi(t) 

X3(t) 

= 

X(10) 

X(13) 

X(16) 

x;(t) 

xjct) 

X3(t) 

X(ll) 

X(14) 

X(17) 

(34) 


Units  of  position,  velocity,  and  acceleration  are  meters,  meters/second, 
and  meters/second  , respectively.  A complete  description  of  the  target  paths 
is  given  in  Appendix  B. 

Path  A. 


♦INSFRT' 

X(14) 

y(17)=-43*SIN(’THE/5) 


X(0) 


1000.0 

1000.0 

0.0 

200.0 

0.0 

0.0 

0.0 

0.0 

0.0 

(35) 
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Path  B 


♦■I'JSERT 

Pir  = 7. 1415926'=: 

X (11)  = -2.6344  51  :328»SIN{2*P  IT) 
X(14)=-10C»SIM(PIT) 


X(0) 


4000 . 0 

0.0 

0.0 

-599.577531 

31.83098861 

0.0 

0.0 

0.0 

0.0 

(36) 


Path  C. 


•INSERT  TflRTAS.9 

IF(TI»^F.LT.2.)  X(ll)  = -3:0 
IF(TI'<»E.GE,2.)  X(ll)=40 


X(0) 


5000.0 

1000.0 

1000.0 

-250.0 

0.0 

0.0 

-300.0 

0.0 

0.0 

(37) 
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Coast  Mode  Simulation 


The  coast  mode  was  simulated  by  defining  a variable  flag,  ICOAST,  which 
was  set  to  FALSE  when  the  normal  mode  was  to  be  simulated  and  set  to  TRUE 
when  the  coast  mode  was  to  be  simulated.  The  original  program  defined  two 
additional  variable  flags,  MIROT  and  IROT,  which  were  utilized  to  check  if 
the  simulated  target  was  detected  in  the  tracker  field  of  view  or  the  tracker 
gate,  respectively.  (Ref  2:28-31) 

Through  the  use  of  multiple  IF  statements,  when  the  coast  mode  was 
selected  (ICOAST  = FALSE),  then  MIROT  and  IROT  were  set  to  FALSE,  even  if 
the  teLTget  was  actually  located  in  the  field  of  view  or  gate,  to  simulate 
sensor  failure.  When  these  two  flags  became  FALSE,  the  tracker  was  forced 
into  a search  mode,  increasing  the  size  of  the  field  of  view  in  order  to 
look  for  the  target. 

The  rate  memory  function  of  the  coast  mode  was  simulated  by  setting  all 
the  pointing  error  samples  to  zero  so  that  the  proportional  plus  integral 
tracker  controller  command  was  held  constant  at  a value  equal  to  the  last 
Integrator  output.  Also  during  the  coast  mode  period,  the  Kalman  filter 
update  equations  were  by-passed  by  an  IF  and  GO  TO  statement  to  avoid 
calculations  with  erroneous  data.  The  following  computer  statements  were 
used  to  simulate  the  coast  mode. 

•IMSERT  TI.7 

CO‘1*iON/COAST/I''OAST 

•I'JSERT  TI.37 

IC3AST=FAL5E 

CCC  1‘ISERT  TVAM  AND  TREAP  HERE  CCC 

TVAN=1.0 
TREAP=5.1 

•DELETE  TI.82 

10  IFdIHE.LT.TVAN.OR.TIME.GE  .TREAP)GO  TO  11 

ICOAST=TRUE 

IROT=FALSE 
MIROTsFALSE 
GO  TO  12 

11  IC0AST=FAL3E 

12  IFCTIHIO.EO.OGO  TO  150 
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•DELETE  TI.d7 

30  IF(  I70T.-D.T?U’^.  aNO.ICOftST  .ET.FALSF)  GO  T^  35 
•DELETE  TI.91 
•helETE  TI.92 
•DELETE  TI.IC: 

IF{  MI^OT.cQ.T^  'f.AND.ICOAST.  EQ.FALEE)  GC  ”^0  100 
•DELETE  ’■I.1C2 

60  rF( IDDAST.FO.-'  L9E) MIR0T=T ’UE 
♦DELETE  ti.178 

IF(TOOA''T.‘^D.FLLr'E)  MIROT=T  ^JJE 
•DELETE  TI.133 

IF(IDDAST.FQ.F^,LSE)I«OT  = TR’tE 
•DELETE  TI.21‘* 

I'^dCOaST.EQ.FALEE)  I«OT  = T‘^')E 
♦DELETE  TI.23‘t 

IP(  ICOA3T.  e-i  LGE)  MI^OrrM  IN3  ( ISMDLL,  I FAMP.) 
•DELETE  TI.2ua 

240  IFdCOAST.ED.FALFE)  GO  TO  241 
ERAZS=0 

ERELS=0 
E,’;A7L  = 0 
ERELL=0 

241  WCA7=A9^ {X(76) •Ae^ERAZS) 

•INSERT  AIDTRK.Z 

COMMON/COAsT/irOAST 
•DELETE  AIDTRK.ESS 
313  IF(ICOAST.ED.-l)  GO  TO  130 
CALL  KALMAN 
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Uncertainties 


The  range,  pointing  error,  and  gyro  bias  drift  uncertainties  were 


simulated  as  follows.  (GAUSS  is  the  subroutine  that  calculated  the  three 


sets  of  normal,  zero  mean  measurement  errors  that  were  added  to  the  range 
and  pointing  error  measurements.) 

RaJig®*  '►I'4S.ERT  AnTRK.l4 

COMMOM  /=!ANGE/  RR(3) 

•DELETE  ftIlTRK.8? 

Rf>(3)=RP(2) 

RRf 2) =RR(1) 

CALL  GAUSSdjl,  0.  ,5.) 

RR(1)=R+X(231) 

•DELETE  AIDTRK.123 

ERD=(3»R°(1) -4  RR(2)+RR(3)  )/(2*DT) 
•DELETE  AnTRK.296 

CALL  GAUSS (1,1»  0 .,5.) 

RR(1)=R+X(231) 

•I'ISERT  APTXVA.12 

COMMON  /RANGE/  RR(3, 

DIMENSION  RR(3) 

•INSERT  APTXVA.196 

IF(ATENA.NE.P.DR.TIM10  .£Q.O)  GO  TD  50  0 
RR(3)=RR(2) 

RR(2)=PR(1) 

CALL  GAUSS(1,1*  0 .,5.) 

RR(1)=R+X(231) 

500  CONTINUE 


(Where  the  units  of  range  measurement  are  meters.) 


Pointing  Error. 


•DELETE  TI.170 

163  CALL  GAUSS  (2,1,0  30001  ) 

AZER=1  00  0.0^(E'»A-^i-X{232)  ) 
•DELETE  TI.175 

call  GAUSS(3,1,  0 .00301  ) 

ELER  = 1GJ3.0*(E'>PL*X(2  33)  )•  TESOS 


(Where  the  units  of  angle  measurement  are  radians.) 


Gyro  Bias  Drift. 

•DELETE  APTOYA.71 

0Y(  l)=WCA7L-»(I“!-AZGC*’4.b48  1368 E-7 
•DELETE  APTOVE.TO 

OY(31>  =WCELL-WIF-ELGC*4.84  81 365 E-7 


(Where  0.1  degrees/hour  = 4,8481368  x 10  radians/second.) 
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Additional  Output 


The  following  statements  were  inserted  in  order  to  provide  additional 


output  data  not  provided  by  the  original  program.  They  are  listed  here 


for  information  only. 

•I'JSERT  riMRSP.73 

6J5  rORMAT {iy,F5.2,5  (£12*4) ) 

505  (IX,  5 (r  1?  . -)  ) 

NUM^FRrl 

JCHEC<=0 

WRITE  {15,&05)T'MF,X(11),X(14),X(17),A'1,F1 
WRITr(  15,^0  6)  E^,'^3,EM,THET  A0,X  (4)  ,X  (S) 

»ISSERT  TI^IRFD.l^l 

IF(Tr’iF.LT.X(81  ) ) GO  TO  ^ 01 
IF( JOHECK.EQ.l) GO  TO  6C1 
SA7=SIN(X(l) ) 

CAZ=COS(X(l)) 

SEL=SIN{X(5) ) 

CEL=COS(X(5) ) 

C11=CAZ^CEL 

C12=-CAZ 

C13=CAZ*SFL 

C21=SA7*CEL 

C22=CAZ 

C23=SAZ*SE. 

C31=-SEL 

C33=CEL 

JC4EC<=1 

601  E1=C11*X(252) ♦^12*X (235) +: 13*X (258) 

E2=C21*X(252)  4-0  2 2*X(2  55)  +C  230*X  (258) 

E3=C31*X(2F2) ♦033*X(258) 

AM=SQRT(X(11)*  2 + X(14)*»  2+  X(17)**2) 

EM=SQRT(F1**2  + E2*»2  + E3*»  2) 

A‘100TEM  = X ( 11)  ♦'"H-X  (14)»E2*-  X(17)  ^FZ 

IF(AOS (AM) .LT.l. OE-IO.OR.A  ^S(EM) .LT.l. OE-19)GO  TO  603 
CTMETA  = Ai-tr)OTEM/  (AM^EM) 

THETA=ACOS(OTHrTA) 

THiTAO=(THETA»130)/3.14l59 26536 
GO  TO  60^- 
603  THETAO=C.O 

6 04  WRITE  ( 15,605)  TIMF,X(  11)  , X(  14)  , X (17)  , AM  ,=•! 

WRITE( 15,606) E?,F3,EM,THETA0,X(4),X (3) 

NUM5ER=NUM5ER4-1 
•INSERT  TI4RSP.215 
607  F0RMAT(1X,I6) 

WRITE(14,60'')  N'lMRER 
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Program  Corrections 


The  following  statements  corrected  programming  errors  that  were 
present  in  the  original  program  available  at  Wright -Patterson  Air  Force 
Base,  Ohio.  They  are  listed  here  for  information  only. 


'►OELETF  ri-tifjo.lft? 

IP(MAT.LE.O)^-n  ■'0  365 
•INSERT  riMINP.197 
385  CONTINUE 
♦INSERT  TIMIMP.R34 
GO  TO  hEO 
♦INSERT  TINHOL.S 
NS=X(208) 

♦INSERT  TINN0L.9 
CALL  APTINI 
♦INSERT  TINMOL.IE 
CALL  APTOYA 
CALL  APTOYE 
♦INSERT  TI-riOL.15 
CALL  APTXVA 
CALL  STATS 
♦INSERT  riNMnL.13 

CALL  APTLIN(TS'P, 0) 

♦INSERT  TIMMDL.21 
CALL  APTNON 
♦DELETE  APTOLK.49 

DATA  I VO /16, 24, 25,26,27,28  .29,3  0 ,132^0  / 
♦DELETE  KALMAN. 23 
♦DELETE  KALMAN. 42 
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Appendix  B 


Simulated  Target  Paths 

Three  representative  target  paths  were  devised.  By  comparing  the 
modified  algorithm  system  performance  with  that  of  the  unmodified  system 
as  each  path  was  tracked,  any  Improvement  or  degradation  due  to  the  modi- 
fication could  be  found. 

The  tracker  base  was  simulated  to  be  stationary  during  all  computer 
runs.  All  target  paths  were  defined  In  relation  to  an  arbitrary  inertial 
coordinate  system  with  the  three  primary  directions  of  Xj^,  and  X^. 

The  origin  of  the  coordinate  system  was  defined  to  coincide  with  the  tracker. 

Path  A 

The  first  simulated  path  was  designated  Path  A,  and  it  represented 
a circular  target  path  with  a 1000  meter  radius  contained  in  a plane 
perpendicular  to  the  X^  inertial  axis  at  a distance  of  1000  meters  from 
the  tracker,  centered  on  the  Xj^  aixls.  The  speed  of  the  target  was  main- 
tained at  200  meters/second  and  the  target  acceleration  was  maintained  at 
a constant  magnitude  of  40  meters/second  . The  target  acceleration  was 
always  directed  toward  the  center  of  the  circle.  The  acceleration  applied 
to  the  target  in  Path  A was  equal  to  approximately  4.08  g. 

Path  B 

The  second  simulated  path  was  designated  Path  B,  and  it  represented  a 
constant  speed,  sinusoidal  target  path  completely  contained  in  the  Xj^X2 
inertial  plane.  At  t = 0,  the  target  was  located  4000  meters  from  the 
tracker  on  the  X^^  axis.  Subsequent  motion  oscillated  at  a frequency  of 
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0.5  Hertz  about  the  axis  toward  the  tracker.  The  speed  of  the  target 
was  maintained  constant  at  approximately  600.42  meters/second.  The 
acceleration  of  the  target  in  the  direction  varied  fx'om  approximately 
-0.27  g to  +0.27  g*  The  acceleration  of  the  target  in  the  X2  direction 
varied  from  approximately  +4.08  to  -4.08  g. 

Path  C 

The  final  simulated  path  was  designated  Path  G.  It  represented  the 
path  of  a target  that  was  accelerated  for  2.0  seconds  along  a straight 
line  by  an  approximate  30.61  g acceleration,  then  decelerated  in  the  same 
direction  by  an  approximate  -4.08  g acceleration.  The  target  was  initially 
located  at  a point  5OOO  meters  in  the  X^  direction  and  1000  meters  in  the 
Xg  and  X^  directions  from  the  tracker.  The  target  moved  in  the  negative 
Xj^  direction,  parallel  to  the  X^  axis. 

Fig  6 illustrates  the  three  target  paths  simulated  during  the  investi- 
gation of  this  thesis.  Initial  conditions  for  each  path  are  listed  in 
Appendix  A under  the  subsection.  Path  Simulation. 
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